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ABSTRACT 


Th« respons* of solar models to perturbations of the 
efficiency of convective energy transport is studied for a 
number of cases. Such perturbations primarily effect the 
shallow superadiabatic layer of the convective envelope (at 
depths i 10^ km below the photosphere) . Independent of the 
details of the perturbation scheme, the resulting change in 
the solar radius (AR/R) is always very small compared to the 
change in luminosity (AL/L) . This appears to 1 true for any 
physical mechanism of solar variability which operates in the 
outer layers of the convection zone. 

Changes of the solar radius have been inferred by Dunham 
et al. (1980,1981) from historical observations of solar eclipses 
in 1715 and 1925. Considering the constraints on concurrent 
luminosity changes, this type of solar variability must be 
indicative uT changes in the solar structure at substantial 
depth'i below the superadiabatic layer of the convective 
envelope. 


Subject headings ; convection - Sun: general - Sun: interior 
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1. INTRODUCTION 


The Sun is often cited aa a paradigm of constancy » both 
in common folklore and in a stricter astrophyiiical context. The 

theory of stellar structure and evolution supports this view by 

placing the Sun in the main sequence stage, where conversion of 

hydrogen to helium provides a stable energy source for some 10^^ 

years. However, there are a number of physical mechanisms which 
can produce short term variations in the Sun that are not 
addressed within the framework of classical stellar evolution 
theory. Some of these mechanisms are: nonradial instabilities 

associated with the distribution of He^ in the solar core 
(Dilke and Go^gh 1972) , beat interactions between g-mode 
oscillations in the core (Wolff 1976) , effects of dynamo- 
generated magnetic fields on the convective envelope (Thomas 
1979; Spiege' and Weiss 1980), and fluctuations in energy 
transport through the convective envelope d”e to the stochastic 
nature of turbulent convection (Dearborn and Newman 1978) . In 
most cases, these mechanisms are too complex to be modeled with 
any confidence on a pur'^ly theoretical basis. For thi^ reason, 
it is advisable to look to the observations for guidance and to 
use ^ hoc modeling to relate observations to batsic physical 
mechanisms . 

The observations we address in this paper concern varia- 
tions of the solar luminosity and radius. We will use perturbed 
solar models to investigate whether these observations can be 
related to fluctuations in energy transport in the convective 
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•nv«lop« of thtt &un. Th« poxtinont observation! are aummarized 
in I XI; the model a are deacribed in fill. The concluaion derived 
in fill ia that aolar model a do not reapond to changea in con- 
vective energy transport in the manner auggeated by the obaer- 
vationa. This auggeata either that stellar evolution models 
are inappropriate for investigating solar variability or that 
som> other physical mechanism is responsible for the variability. 
These alternatives are considered in some detail in iSIV and V, 
respectively. 
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II. A SUMMARY OF THE OBSERVATIONS 
a) Irradlance Measuremanta 

Meaaurei^ients of the solar luminosity actually refer to the 
irradiance S (solar constant ) , which is the flux emitted in the 
direction of the Earth. Ground-based measurements were made by 
the Smithsonian Astrophysical Observatory from 1902 to 1962, 
with the mohrt extensive data covering the period 1923 to 1954 . 
Based on a statistical analysis of the monthly means for the 
latter period, Hoyt (1979) concluded that the measured varia- 
tions (of order +0.5%) are due primarily to incomplete compen- 
sation for atmospheric effects (see, also, Roosen, Angione, 
and Klemc)ce 1973) . According to Hoyt, the Smithsonian data 
can probably set a limit | A S/S | ^ 5 x 10 for irradiance 
variations over the period covered by the measurements. 

Because of the large corrections for atmospheric absorp- 
tion (typically 25 to 40%) required in ground-based observations, 
modern measurements have been made at balloon altitudes or higher. 
Frohlich (1977) and Frohlich and Brusa (1981) have reviewed the 
balloon and rocket measurements and conclude that |AS/S| < 1.5 x 
lO”^ from 1969 to 1980 (i.e., over one sunspot cycle). The sub- 
ject of solar irrfiidiance variations is often controversial and 
the claimed detections of variability (e.g., Kondratyev and 
Nikolsky 1970) must be treated with healthy skepticism. 

The Nimbus 7 and Solar Maximum Mission (Sf4M) satellites 
have recently detected variations in the solar irradiance 
(Hickey et al . 1980; Willson et al . 1981) . The reality of these 
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variations is confirmed by their simultaneous detection by the 

two independent experiments. For the two-year interval (1978-80) 

covered by the available date, the irradiance varied by as much 

as 0.2% from the mean (|()S/S| • 2 x 10~^) . These fluctuations 

are rapid (time scales < 10 days) and upper limits on secular 

trends are much smaller. Perhaps more important, Willson et al . 

(1981) found that the major variations in the SMN measurements 

coincided with the passage of large nunspot groups across the 

central meridion of the solar dislc. This suggests that the 

observed irradiance variations are due to local bloc)cing of t)>e 

solar radiation in the vicinity of active regions. In this 

case, caution must be used in relating observed irradiance 

variations to changes in the solar luminosity. A recent 

analysis of the SMM data by Oster, Schatten, and Sofia (1981) 

concludes that the irradiance variations are dominated by 

anisotropy effects (directionality of the radiation) and that 

true variations of the solar luminosity are smaller than one- 

tenth of the irradiance variations (i.e., AL/L ^ 0.1 (AS/S) 

-4 

^ 2 X 10 , for short-term fluctuations) . In this case, 

observations over a long period of time will be required to detect 
global luminosity va iations. 

b) Radius Measurements 

Solar radius changes of sufficient amplitude would show up 
in the transit measurements made at a number of national obser- 
vatories. The measurements made at the Greenwich Royal Obser- 
vatory between 1850 and 1937 have been analysed by Sofia et al . 
(1979). They found that jAR/R| < 2.5 x 10 ^ for any secular 
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trend during this period.^ Shapiro (1980) used timings of the 


'^Eddy and Boornazian (1979) also used the Greenwhich data 
to loo)c for trends in the solar radius and claimed a steady 
decrease amounting to /^R/R ■ -1 x lO”^ for the period 1836-1953. 
However, this large rate of shrin)(age is ruled out by a number 
of other measurements (described in the text) . 

passage of Mercury across the solar dis)f to establish an indepen- 
dent limit on long-term changes in the solar radius. For the 
period 1736-1973, Shapiro found |AR/R| < 3.6 x lO"^ for any 
linear trend. It should noted, however, that both the 
transit data and the Mercury data show large amounts of scatter 
so nonsecular variations of the radius could be much larger 
than the limits given above. Gilliland (1981) has reanalyzed 
transit measurements and the Mercury data. By combining several 

data sets, he found marginal (2o) evidence for a 76-year oer- 

-4 

lodicity, with an amplitude AR/R ^ 2 x 10 

The solar radius can be very accurately determined from 
observations of the width of the path of totality in a solar 
eclipse. Unfortunately, the required data (eclipse duration 
near the edge of the path of totality) are only available for a 
few eclipses. The inferred radius changes, relative to a stan- 
dard reference value, are listed in Table 1. It should be noted 
that the uncertair in the 1715 value refers to an estimate of 
the maximum possible error since insufficient data are avail- 
able to form a statistical error estimate (which would be 

significantly smaller) . The 1715 and 1925 eclipses give radius 

-4 

values which differ by AR/R ~ 6 x 10 from recent values. The 
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diffcr«nc« between th& 1925 and 1980 valuas it highly signifi- 
cant (graatar than 5o) . 

Conaidaring tha datas of tha obsarvations, it ia apparant 
that thara is no linaar trand, and tha changas may wall ba 
stochastic. Tha 1715, 1925, and 1976 measuraments all occur 
on tha rising branch of tha sunspot cycla so thara is no 
obvious correlation with solar activity. Also, according to 
the 76-year periodicity claimed by Gilliland (1981), the 1925 
and 1980 eclipses should have occurred near maximum solar 
radius and the 1715 eclipse should have fallen near minimum 
radius. This is clearly in conflict with the data, indicating 
that Gilliland's suggested periodicity is not real. The ques- 
tion of periodicity, and time -dependent behavior in general, can 
only be addressed when additional observations become available. 

At this point, it is only possible to determine that radius 

-4 

changes AR/R ^ x 10 occur on some time scale t < 50 years. 

c) Photospheric Velocity Fields 

Solar radius changes AR, on some time scale t, imply expan- 
sion velocities V AR/t, which would cause a variable center- 
to-limb Doppler shift across the solar disk. If we assume that 
the Sun is presently in a variable state, then modern Doppler 
measurements can be used to set upper limits on AR/t. 

Figure 1 shows the resulting limits on the time scale of 
solar radius variations. The limit of spectroscopic detection 
has been set at 1 m/s, on the basis that large-scale variable 
velocity fields as small as 2 to 3 m/s have been measured 
(Severny, Kotov, and Tsap 1976; Howard and LaBonte 1980) . This 
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detection limit eete a lower limit of -B days for the time 

-4 

scale of recent radius changes of order AR/R ■ 6 x 10 . Thus, 

reasonable (though not certain) limits on the time scales of 
the radius changes are: 5 days < t < 50 years. While not very 

useful from an observational viewpoint, these limits are impor- 
tant in defining the modeling regime (see 5§IZI and IV). 

Figure 1 also shows the approximate locations in the V-t 
plane of the 5-minute and 160-minute solar oscillations. These 
oscillations have amplitudes which are factors of 60 or more 
smaller than the radius changes we are considering in this 
investigation. 
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in. PERTURBED SOLAR MODELS 

Ulrich (1975) introduced the notion of modeling eoler 
variability by changing the characteristic length scale in the 
mixing length formalism fer convection. Since this length scale 
is usually specified by a, the (constant) ratio of mixing length 
to the local pressure scale height, we will refer to this pro- 
cedure as an a-perturbation. If an a-perturbation is intro- 
duced in a time-series of stellar evolution models, the re- 
sponse of the models may mimic the response of the Sun to 
thermal fluctuations within the convective envelope. As 
pointed out by Dearborn and Newman (1978), changing the 
mixing length is equivalent to changing the efficiency of con- 
vective heat transport. 

The response of a solar model to an a-perturbation has been 
described in detail by Dearborn an<* Blake (1980) and Sweigart 
(1981). Briefly, for a positive perturbation (increasing a), the 
efficiency of convection increases, causing a sudden increase in 
both L and R. The luminosity perturbation decays slowly, on a 
time scale of order 10^ years. By contrast, the radius pertur- 
bation appears to decay much more rapidly. This is because the 
initial expansion is restricted to the outer, superadiabatic 
layers, while the deeper layers are slowly contracting. Because 
the outer layers have a short thermal relaxation time, the radius 
returns to its unperturbed value after a few hundred years. 
Thereafter, the slow contraction of the inner layers dominates 
and the photospheric radius is smaller than the unperturbed 


value. 
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On* probl«in with such perturbation procedure ia that it ia 
difficult to know what conatitutea a realiatic amplitude for 
changea in a. Aa a reault, change* R and L of almoat any 
aiae can be produced. Sofia ^ (1979) auggeated character- 

izing a particular type of perturbation by 

W - , (1) 

LL/L 

From the behavior of a aolar model deacribed above, it ia clear 
that W will be time-dependent. However, during the firat 50 
yeara or ao after a perturbation, the value of W doea not 
change much from ita initial value. More important, the value 
of W ia, to firat order, independent of the perturbation ampli- 
tude. 

Baaed on model calculations, Sofia ^ (1979) eatimated 

W - 7.5 X 10“^ for an a-perturbat ion. Thia retult was questioned 

by Dearborn and Blake (1980) and Gilliland (1980) , who obtained 

much smaller values for W. A recalculation using essentially, 

the same stellar evolution code as used by Sofia et a^. showed 

that the earlier calculation was, in fact, in error (Twigg and 

Endal 1981) . Values of W resulting from o-perturbations have been 

calculated by Dearborn and Blake (1980) , Gilliland (1980) , Sweigart 

(1981), and Twigg and Endal (1981), using independent stellar 

evolution codes. The reported values (immediately after a per- 

-3 -4 

turbation) range from 6 x 10 to 6 x 10 . Although the range 

of these values shows that the results are fairly code-dependent, 
a general conclusion that W << 1, for an Q-perturbation, is war- 
ranted. We can compare this result with the value of W implied 
by the observations described in III. 
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Th« •clips* m«ssurttn«nts of thm solsr radius (fxi.b) show 
that solar radius changas of ordsr dR/R ■ € x 10*^ hav« oeeurrod 
witVin this csntury. On tha othar hand* tha programs tot datac- 
ting solf^r irradianca changas (fXZ.a) hava failad to dataet any 
changas. Thasa programs covar tima scalas ranging frcm 60 yaars 
dov^ to a fraction of a day. Tha obsarvations raquira Us/S| < 

5 X 10*^. Baaad on tha axpactatior that |0L/L| |08/8|» this 

givas 

|W| > 0.1 (2) 

Thus, tha obsarvations raquira that tha magnituda of H ba 
significantly graatar than tha valuas obtainad from tha 
thaoratical calculations. This suggests that aithart (a) 
although tha observed radius changes may be due to chances in 
efficiency of convection, the real variability cannot be modeled 
as a simple linear perturbation; (b) the approximations made in 
stellar evolution calculations (hydrostatic equilibrium, mixing 
length theory of convection, etc.) are not appropriate for this 
type of modeling; or (c) the observed radius variations are not 
associated with changes in the efficiency of convection in the 
solar envelope. The first possibility can be tested by using 
models to estimate the importance of nonlinear behavior, as 
described below. Discussions of (b) and (c) will be deferred 
to iilV and V, respectively. 

If the solar radius and luminosity respond linearly to an a- 
perturbation, W will be independent of the amplitude of the per- 
turbation. In this case, the value of W will not depend on the 
amplitude or on the time interval between successive perturbations 
(Since the cumulative changes AR and AL add linearly). However, 
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th« thermal atxuctura of the Sun la govarnad by nonlinear aqua- 
tlona and la characterized by a broad range of relaxation tlmea. 
Under theae condltlona, It la not obvloua that all poaalble 
valuea of W are obtained by the alngle, email -amplitude pertur- 
batlona uaed In prevloua calculatlona. More complex perturba- 
tion achemea (perhaps Involving a-perturl itiona of models which 
are not in thermal equilibrium) might give rise to lai'qer valuea 
of |W|, more in accord with the observationa . Such perturbation 
achemea may alao give a better repreaentation of the real 
aituation, in which the Sun ia never in true thermal equilibrium. 

To teat thia poaaibility* we calculated the reaponae of a 
solar model using various schemes for perturbing the mixing 
length in the convective envelope. The calculations were 

9 

started from a 1 model evolved to an age of 4.7 x 10 yr, 

with the input physics described by Endal and Sofia (1981) . 

To minimize numerical noise* a large number (700) of interior 

zones was used and the convergence criteria for all internal 

iterations were tightened by several orders of magnitude beyond 

normal stellar evolution values. In addition* the difference 

equations were rewritten in centered ("conservative”) form to 

increase numerical stability at short time steps. The outer 

boundary conditions were provided by a fine grid of static 

•■ft 

envelopes of mass 6 x 10 The thermal relaxation time (1.4 

years) of this mass effectively limited the time resolution of the 
calculations. 

One way to test for nonlinear effects is to test for ampli- 
tude-dependence of W. Figure 2 shows the value ot W for pertur- 
bations of various amplitudes imposed on an equilibrium model. 

The solid line refers to values obtained using a standard miring 
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length formalism, with the mixing length proportional to the 
pressv«‘e scale height. The dashed line refers to an alternative 
formulation, idiich will be discussed in ilV.d. In both cases, 

H varies by « 101 over a factor of 40 range in perturbation 
amplitude, indicating nearly linear behavior. This is consistent 
with the linear relationship between Aa and AL found by Dearborn 
and Newman (1978) . The smooth behavior in Figure 2 indicates 
that numerical noise was under control in the calculations. 

We can also test for nonlinear behavior by imposing a 
series of perturbations on an evolutionary sequence of models 
(with short time steps) . If the time intervals between succes- 
sive perturbations are shorter than the relaxation time of some 
part of the model, then the perturbations may add in a nonline r 
fashion since they are applied to models which are not in thermal 
equilibrium. Although this effect could be examined using a 
purely random or arbitrary sequence of curbations, it seems 

more instructive to concentrate on some Icnown feature of solar 
variability, such as the 22-year magnetic cycle. For this pur- 
pose, a 1%-amplitude, 22-year cycle of a-pertur!oations was imposed 

2 

on a sequence of models spaced at one-year intervals. The 
evolution was followed for 374 years (17 cycles) . In order to 


2 

Using time steps shorter than the 1.4-year resolution 
allowed by the mass of the static envelope means that the time 
domain was over-sampled. This does not create a problem since 
the time scale for the a-perturbations was much longer (22 years) . 


simulate the noisy character of the real solar cycle, the pertur- 
bations were generated by passing a sequence of random numbers 
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through a narrow bandpass filter centered at 1/22 cycle per year 
(see Barnes, Sargent, and Tryon 1980) . 

The percentage a*perturb^. cions 


Ma^) 



X 100 


(i) 


are shown in Figure 3. Here, is the mixing length to pressure 

scale height ratio used in mod' i (i « 1, 374) and gives 

the unperturbed value. Applying the perturbations to the models 
yielded the following quantities of interest; 


«(R^) 



X 100 , 


(4) 


and 


<* (L. ) 



X 100 , 


W. « (6R.)/(6L.) . 


(5) 


( 6 ) 


These quantities are shown in Figure 4. From Fioure 4(c), it 
is apparent that the values of for multiple perturbations do 
not differ greatly from the value (W * 6 x 10”^) produced in our 
models by a single perturbation.^ 


^The curves in Figures 4(c) and 5(c) show occasional sharp 
spikes. These spikes result from slight phase shifts of the 6(R^) 
and 6(L^) curves. Considering the limited time resolution of our 
calculations, these phase shifts cannot be considered real. On 
the other hand, the broad dips in appear to be real. This dips 
invariably occur in the negative parts of the ct-cycle and may be 
associated with the thermal relaxation process. 
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j^t this point, w« digrosB ilightly to consldor an obsorvational 
problain: tha observations yield c!>anges AR - where neither 

R^ nor Rj can be identified with the equilibrium radius used as a refer*- 
ence in the theoretical models. To examine the effect of using 
different reference values to define we replaced the equili- 
brium values R^ and in equations (4) and (5) by values from 
n previous time steps, i.e.. 



where x refers to R and L and equation (6) still defines 
An example of the effect of altering the reference value in 
this ...diiner, for n«ll, is shown is Figure 5. These values were 
computed from the same runs as shown in Figures 3 and 4. 

Although the 6(R^) and 6(L^) curves are now considerably more 
noisy, the general level of the values is not substantially 
altered. Curves for n > 1,2,5, and 22 were also examined, giving 
similar results. (The n ■ 11 case shows the largest deviation 
from Figure 4.) The stability of the W values for a wide variety 
of perturbation schemes and operational definitions indicates 
that this parameter does, indeed, provide an excellent character- 
ization of the a-perturbation. 

Finally, it is possible that the a-perturbations imposed on 
our models were too rapid to allow nonlinear effects associated 
with the relaxation process to arise in the deeper parts of the 
convection zone. To test this, a smaller number of cycles was 
computed with periods of 44 and 110 years. The results confirm 
the conclusions given abo.e. Time scales longer than 110 years 
are of doubtful relevance to the available observations. 
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We conclude that a-perturbationa in stellar evolution models 
of the Sun will not yield values of W in the allowed observational 
range, irregardless of how these perturbations are introduced . The 
basic problem is that it is not possible to produce significant 
radius changes without introducing unreasonably large variations 
in the luminosity. For W » 6 x 10~^, the observed radius changes 
AR/R ^ 6 X 10 imply a lOt variation in luminosity. Aside from 
the direct observational constraints, such a large luminosity 
variation would produce easily detectable changes in the earth's 
climate (cf . Budyko 1969) . For the even smaller values cf W 
found by Gilliland (1980) and Sweigart (1981) , the problem is 
more severe. This suggests that either t*’? response of the Sun 
to changes in convective efficiency is incorrectly modeled in 
the calculations described above, or that some other physical 
mechanism is responsible for the observed radius changes. 
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IV. ANALYSIS OF THE MODEL APPROXIMATIONS 

Thtt approximations made in stellar evolution calculations 
have been thoroughly described in a number of standard texts 
(e.g., Schwarzschild 1958; Cox and Guili 1968). However, calcu- 
lations such aj those presented in 6 III are very different from 
normal stellar evolution calculations. For this reason, the 
important approximations, and their justifications in the present 
context, should be examined in some detail. 

a) Spherical Symmetry 

The full amplitude of the observed radius changes. AR/R 
6 X lO”^ is more than 10 times larger than the distortion of the 
solar disk due to the permanent oblateness. " **pole^^** * 

5 X lO”^ (Dicke and Goldenberg 1967; Hill and Stebbins 1975) . and 
30 times greater than the periodic distortion. (R^^ - - 

2 X 10“^. reported by Dicke (1976) . Unless the oblateness measure 
ments (in 196G and 1973) refer to years of anomalously small dis- 
tortions. we can conclude that departures from spherical symmetry 
do not play a role in the observed radius variations. 

b) Quasi-hydrostatic Equilibrium 

The usual statement in stellar evolution studies is that 
hydrostatic equilibrium is valid if the relevant time scales are 
much longer than the free-fall time: 

— 5“"— • 

■ 2/ R /GM 2 . 50 minutes, for the Sun. 


( 8 ) 
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In the present context, we can quantify this statement. Any 
radius change must be due to an inbalance of the gravitational 
and pressure gradient forces. Let f be the fraction of the pres- 
sure gradient force not balanced by local gravity. Then, assum- 
ing spherical symmetry, 




fg » 


(9) 


where g is the gravitational acceleration. For the small radius 
changes considered here (AR << R) , we can assume that g is essen- 
tially constant and let f denote a mean value over an expansion/ 
contraction phase. Successive integrations of equation (9) give 


,, dR 
V s ——.I 

dt 


fg t 


( 10 ) 


and 


A5 , 1 is t2 

R 2 R 



( 11 ) 


where t is now the duration of an expansion/contraction phase. 
Fzom equation (11) , we obtain 


f 



( 12 ) 


In principle, f can be very small and still violate the 
hydrostatic assumption, since strict hydrostatic equilibrium 
implies f * 0. Quasi-hydrostatic equilibrium includes the rates 
of change of gravitational potential energy (U^) and internal energy, 
but ignores the rate of change of )cinetic energy (Uj^) . This appro- 

• s 

ximation is valid if Ug >> Uj^. Using equations (9) to (11) to 
evaluate these rates for a unit mass at the solar surface gives 
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With equation (12) and Z)R/R ^ 
becomes 


t (13) 

f^g^ t . (14) 

-4 

, the ratio of these rates 



For any nearly homologous change, this ratio is smallest near the 
surface so the quasi-hydrostatic approximation is amply justified 
for t > t^^. If the lower limit on the time scale derived in §II.C 
is applied, then t > 150 t^^ and Uq/Uj^ > 7 x 10^. 


c) Approximations for Radiative Transfer 
In a stellar interior, the diffusion approximation gives an 
excellent representation of energy transport by radiation 
(Schwarzschild 1958, Chapter 2) , so this is not normally considered 
a significant source of error. However, the present investigation 
is specifically concerned with determining the radius of the Sun, at 
the photosphere. In this region, the diffusion approximation 
brealcs down and it is necessary to use a different scheme for 
integrating the structure and locating the photosphere. The cal- 
culations described in §III used the Eddington approximation (see 
Paczynslci 1969) for optical depths t < 2/3, and the diffusion 
approximation for t ^ 2/3. We can chec)c this scheme by comparing 
the atmospheric structure in our equilibrium model with an empirical 
solar atmosphere, derived from the observed solar spectrum. Figure 
6 shows such a comparison with the empirical solar atmosphere 
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derived by Vernazza, Avratt, and Loaaar (1976) . The quantitiea 

shown are: T(t) - tsinpsrature vs. optical dspth, and X(t) - 

physical depth (with respect to t ■ 2/3) vs. optical depth. 

The latter quantity determines the radius assigned to our 

models. Our X(t) deviates from the empirical relationship by 

less than 30 )(m over the primary atmospheric region 0.01 ^ t ^ 

10. This should he compared to the radius changes being modeled — 

-4 

dR ^ 6 X 10 R B 418 )(m. Thus, the Eddington approximation is 
sufficiently accurate for our purposes. 

d ) Mi xing Length Theory of Convection 

The mixing length theory of convection has been vigorously 
criticized in the literature. In particular, Gilliland (1980) 
has argued that uncertainties in the treatment of convection 
invalidate present attempts to model solar variability. While 
such criticisms are easily made, the fact remains that a mixing 
length formalism provides the only presently feasible Means of 
carrying out such calculations. In this case, an assessment of 
the uncertainties introduced by this approximation is useful. 

Because there is no "absolute" theory to compare with, this 
assessment will necessarily be qualitative. It is conv'*nient 
to consider first the uncertainties in the equilibrium structure 
and then the uncertainties due to time-dependent effects. 

Any reasonable convection theory will predict that the tempera- 
ture gradient in the bulk of the convection zone is nearly adiabatic. 
As emphasized by Gough and Weiss (1976) , the adiabat is determined 
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by the requirements that the base of the envelope join smoothly 
to a radiative interior with a nuclear luminosity of 1 and 
that the radius of the model be 1 Our equilibrium models 

match these parameters to better than 1 part in 10^, and this 
fixes the mean adiabat with very little dependence on the con- 
vection theory. 

In a local mixing length theory » significant deviations from 
the adiabatic gradient are found in a shallow transition layer, 
with a depth of ~10^ km. The superadiabatic structure of this 
layer is determined by radiative losses from the convectinn 
elements. Since different mixing length theories assume different 
bubble geometries, the radiative losses and the structure of the 
transition layer depend on the specific form of the mixing length 
theOi.y being used. These specific forms are more or less arbitrary, 
so this introduces an uncertainty in the equilibrium structure 
(near the transition layer) . As shown by Dearborn and Blake (1980) , 
the direct effect on energy transport due to an a-perturbation is 
restricted to the transition layer so it would appear that uncer- 
tainties in the stru 'ture of this layer could translate into a 
major uncertainty in calculating W for an a-perturbation. 

To estimate this uncertainty, we calculated the response of 
a solar model to an a-perturbation using two fairly different 
mixing length theories. As a standard theory, we used the formu- 
lation of B&hm-Vitense (1958), as modified by Henyey, Vardya, 
and Bodenheimer (1965) . This formulation assumes that n bubble 
moves a distance (the mixing length) equal to its own diameter. 

The mixing length is assumed to be proportional to the local 
pressure scale height. As an alternative model, we used the 
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formilation dascribad by Chan, Wolff, and Sofia (1981). In thia 
caaa, tha mixing langth is chosan to ba proportional to tha dan- 
aity scala haight in tha bubbla and tha bubble geometry is 
chosan to maximize tha convective flux. Both theories contain 
an arbitrary length scale which must be calibrated by matchincr 
the solar radius and luminosity, as described earlier. Figure 
7 shows the predicted departures from the adiabatic gradient in 
the solar transition layer. Note that the degree of super- 
adiabaticity in the two models is very different, though the 
depth of the transition layers are nearly identical. The latter 
agreement is a direct result of choosing mixing lengths to match 
the observed solar radius. The responses of these models to a 
single a-perturbation were shown earlier, in Figure 2. vhe cal- 
culated values of W differ by only ~20%. This close agreement 
is a direct result of the calibration procedure, which fixes 
the depth of the transition layer. We conclude that, because of 
the observational constraints on the structure of the convection 
zone, uncertainties in the details of local mixing length theories 
do not seriously affect the calculated W. 

The basic reason that W is so small for an a-perturbation is 
that the radius change results from expansion cf the transition 
layer and atmosphere and these regions are very shallow. Nonlocal 
mixing length theories may require a superadiabatic gradient in 
a boundary layer at the bottom of the convection zone. An expan- 
sion of this layer would lead to a large radius increase since the 
entire convection zone would be lifted and expanded (due to the 
increase in the hydrostatic pressure scale height) . However, any 
departure from the adiabatic gradient at the bottom of the convec- 


24 


tion son* ifould d«er*ss« ths total dapth of tha eonvactlon sona. 
Local mixing langth thaoriaa pradict a total dapth of *165,000 km. 
This ia slightly laas than tha minimum dapth of 175,000 km obtainad 
by Rhodaa, Ulrich, and Simon (1977) from an analysis of tha solar 
5-minuta oscillations. A subatantial svparadiabatic gradient 
at tha bottom of tha convection son* would increase this dis- 
crepancy. 

Finally, non* of tha convection theories considered above 
include any dynamical effects. Such affects inay altar tha 
response of the modal to a perturbation if the perturbation 
occurs on a time scale shorter than or comparable to tha con- 
vective turnover time. The observed turnover times at the 
photosphere vary from *7 minutes for granulation to *1 day for 
supergranulation. Deep in the convective envelope, turnover 
times may be on the order of 1 month, but we have already noted 
that this region is nearly adiabatic and difficult to perturb 
by changes in energy transport. While it is not possible to rule 
out nonlocal or time-dependent (dynamical) effects, it appears 
unlikely that such effects would lead to the factor of ten or 
more increase in |w| required to bring the models into agreement 
with the observations. 


e) Code Dependence 

The published values of W for a-perturbations vary by roughly 
a factor of 10, which is comparable to the discrepancy between 
the models and the observations. However, with the exception 
of the incorrect value obtained by Sofia ^ (1979) , all of 

the predicted values of |w| are much smaller than the allowed 
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ob»«rvational rang#. Numarical axparimanta show that ths valua 
of |W| is quits ssnsitivs to ths input physics (opacitiss, for 
instancs) and dstails of ths codss such as interpolation pro- 
csdurss. Howsvsr, while many of these changes produce values of 
W smaller than those shown in Figure 2 , nc reasonable changes 
produce values substantially larger. This is because, as the 
predicted radius change for a given a-perturbation becomes 
larger, it also becomes less sensitive to details of the computino 
method. The predicted luminosity change is very imsensit.i ve to 
details of input physics and coding. 


n 


V. CONCLUSIONS 

Th« tolar radius changaa found from aclipta obaarvationt hava 
tha propartiaa 

AR/R . 6 X 10"^ , |AL/L| < 5 x 10‘^ 

and occur on aoma tlma tcala batwaan 5 days and 50 yaara. Tha 
limit on tha ratio of radius to luminosity changa ia |w| > 0.1. 

mt 

With t! aaa charactariatica, tha radius changaa cannot ba dua to 
fluctuations in energy transport through the shallow auperadiabatic 
layers of the solar convective envelope. Sofia and Chan (1981) 
have examined the response of a solar model to perturbations in 
total pressure due, for instance, to changing magnetic fields. 
Again, they found very small (.10 ) values of W if the pertur- 

bation is confined to shallow layers. These results can be 
understood by noting that a perturbation at some depth in the 
envelope affects primarily the layers at and above that depth. 

The densities and pressures below the perturbed layer are too 
great for these lower layers to be moved any appreciable distance. 
Therefore, any mechanism which is confined to shallow layers will 
produce very small radius changes. We conclude that the observed 
radius changes are large enough that they must be due to mechaniums 
which operate at substantial depths below the transition layer of 
the convective envelope. Sofia and Chan (1981) found W 0.1 for 

pressure perturbations at depths corresponding to (1 - M_/M^) > 

.4 

10 . In our models, this corresponds to physical depths 


X > 27,000 km. 
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At present, the subject of vsristions In globsl solsr struc- 
ture is rich in theoreticsl speculstion but poor in obsorvstionsl 
dsts. Further progress will require more extensive end better 
observations. Simultaneous detection of radius and luminosity 
variations would replace the present lower limit onjwjwith a 
definite value (or range of values) . It may well be true that V' 
is determined by the depth at which the changes originate but 
is largely independent of the specific physical mechanism. In 
this case, additional information (such as time scales, period- 
icities, and possible correlations with solar act vity) will be 
required to pin down the physical mechanisms. Until such obser- 
vations become available, our understanding of the forces which 
affect solar structure on human time scales will lemain vague. 

We thank Drs. K. L. Chan, S. Sofia, and A. V. Sweigart for 
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TABLE 1 


SOLAR ECLIPSE DETERMINATIONS OF THE SOLAR RADIUS^ 


Eclipse 

date 

(AR/R) 

X 10^ 

Reference 


3 May 1715 

+5.4 

+ 

2.l‘' 

Dunham 

et 

al . 

(1980) 

24 Jan. 

1925 

+6.2 

+ 

0.8® 

Dunham 

et 

al . 

(1981) 

23 Oct. 

1976 

-2.4 

+ 

1.5® 

Dunham 

et 

al. 

(1980) 

26 Feb. 

1979 

-0.8 

+ 

0.9® 

Dunham 

et 

al . 

(1980) 

16 Feb. 

1980 

-0.3 

+ 

0.4® 

Dunham 

et 

al. 

(1981) 











0 

The quoted changes are measured with respect to a reference 
radius used in the data analysis. 

b 

maximum uncertainty (not lo) . 

Q 

lo uncertainty. 
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FIGURE CAPTIONS 

PIG.l. -- Th« diagonal line shows the expected expansion/ 
contraction velocities due to radius changes, as a function of the 
time scale t. The intersection of this line with the spectro* 
scopic detection limit (horizontal line) sets a lower limit on 
the time scale of the radius changes. The upper limit is based 
on the interval over which radius changes have been detected. 

The symbols (A) indicate che approximate positions in this 
diagram of solar pulsations with periods of 5 minutes (cf. 

Claverie et 1979) and 160 minutes (cf . Severny ^ 1976) . 

FIG. 2. — The value of W for single a-perturbations, as a 
function of perturbation amplitude. Solid line and heavy dots — 
standard mixing length theory; dashed line and triangles-- 
alternative mixing length theory (see §IV.d) . 

FIG. 3. — Percentage deviation of a from the equilibrium 
value in the 22-year solar cycle simulation. 

FIG. 4. — Changes in the models resulting from the 22- 
year a-perturbations (see Fig. 3) . The curves refer to devia- 
tions of R and L from the equilibrium values. 

FIG. 5. -- Same as Fig. 4, but with 6(R) and 6(L) referring 
to changes in models separated by an 11-year interval. 

FIG. 6. — The T(t) and X(t) relationships for a solar atmos- 
phere based on the Eddington approximation (solid lines) and an 
empirical solar atmosphere (dashed lines) . 
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FIG. 7. -- Difference between the actual temperature gradient 
7 ■ (d&nT/dinP) and the adiabatic gradient 7^^ ■ 0£nT/3'.nP)^^ 
aa a function of depth below the photosphere (t - 2/3) . The 
solid line shows an equilibrium solar model computed with the 
Bohm-Vi tense (1958) mixing length theory and the dashed line 
shows the same model computed with the Chan ^ (1981) 

mixing length theory. 
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